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Abstract 

We investigate the production of gravitinos in a cosmological background. Gravitinos can be pro- 
duced during preheating after inflation due to a combined effect of interactions with an oscillating 
inflaton field and absence of conformal invariance. In order to get insight on conformal proper- 
ties of gravitino we reformulate phenomenological supergravity in SU{2, 2|l)-symmetric way. The 
Planck mass and F- and D-terms appear via the gauge-fixed value of a superfield that we call 
conformon. We find that in general the probability of gravitino production is not suppressed by 
the small gravitational coupling. This may lead to a copious production of gravitinos after infla- 
tion. Efficiency of the new non-thermal mechanism of gravitino production is very sensitive to the 
choice of the underlying theory. This may put strong constraints on certain classes of inflationary 
models. 



^ Onderzoeksdirecteur, FWO, Belgium 



1 Introduction 



The possibility of excessive production of gravitinos is one of the most comphcated problems 
of cosmological models based on supergravity. Such particles decay very late, and lead to 
disastrous cosmological consequences unless the ratio of the number density of gravitinos 
77,3/2 to the entropy density s is extremely small. For example, the ratio of the number 
density of gravitinos to the entropy density s should be smaller than O(10~^^) for 
gravitinos with mass O(IOO) GeV |^. The standard thermal mechanism of gravitino 
production involves scattering of particles at high temperature in the early universe. To 
avoid excessive production of gravitinos one must assume that the reheating temperature 
of the universe after inflation was smaller than 10^ — 10^ GeV [jl], |[|. 

However, gravitinos can be produced not only in the thermal bath after reheating, 
but even earlier, during the oscillations of the inflaton field at the end of inflation. We 
already know that bosons as well as spin 1/2 fermions can be copiously produced by the 
coherently oscillating inflaton field. Quite often this effect occurs in a non-perturbative way 
during the stage of preheating ^. Similar effect may occur for gravitinos. According 
to [P], nonthermal gravitino /mo dull production may rule out certain classes of inflationary 
models which otherwise would be quite legitimate. 

The theory of the cosmological gravitino production is very complicated. Recently 
production of transversal gravitino components (helicity 3/2) was studied in Ref. (see 
also 10 where an attempt has been made to study this question using perturbation the- 
ory). In this paper we will investigate production of all gravitino components, transversal 
(helicity 3/2) and longitudinal (helicity 1/2). As we will show, the production rate of the 
longitudinal gravitino component can be much greater than that of the transversal gravitino 
components.^] 

Since gravitino is a part of the gravitational multiplet, one could expect that their 
production must be strongly suppressed by the small gravitational coupling. Indeed, usually 
the production of particles occurs because their effective mass changes nonadiabatically 
during the oscillations of the inflaton field [Q. This is the main effect responsible for the 
production of the gravitinos with helicity 3/2. The gravitino mass m at small values of the 
inflaton field is proportional to Mp'^W, where is a superpotential. Thus the amplitude 
of the oscillations of the gravitino mass is suppressed by Mp^. That is why production of 
the gravitino components with helicity 3/2 is relatively inefficient P, 0]. 

There is another mechanism to be considered, which is related to breaking of con- 
formal invariance ||^. It is well known, for example, that expansion of the universe does not 
lead to production of massless vector particles and massless fermions of spin 1/2 because 
the theory of such particles is conformally invariant and the Friedmann universe is confor- 
mally flat. Meanwhile, massless scalar particles minimally coupled to gravity (as well as 

"'^In flat spacetime, transversal components, which exist even if supersymmetry is unbroken, correspond 
to the hehcity 3/2, while longitudinal components with non- vanishing iJjq and j^ipi correspond to the helicity 
1/2. Although the helicity concept has a less precise meaning in FRW metrics, we still will use this loose 
definition as a shortcut to the transversal and longitudinal components. 
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gravitons) are created in an expanding universe because the theory describing these parti- 
cles is not conformally invariant. The rate of scalar particle production is determined by 
the Hubble constant H = ^ = Mp^^Jp/3, where p is the energy density. If similar effects 
are possible for gravitinos, they could be much stronger than the effects discussed in the 
previous paragraph. Indeed, H is suppressed only by the first degree of Mp^. As a result 
H typically is much greater than 7713/2 after inflation, so its time-dependence may lead to a 
more efficient particle production. 

The issue of conformal invariance of gravitinos is rather nontrivial, and until now it 
has not been thoroughly examined. For gravitinos with helicity 3/2 the effects proportional 
to H do not appear, and therefore the violation of conformal invariance for such particles 
is very small, being proportional to their mass 0. However, as we will show, the theory of 
gravitinos with helicity 1/2 is not conformally invariant, and therefore such particles will 
be produced during the expansion of the universe even if one neglects their mass. 

But the most surprising effect which we have found is that the production of the 
gravitinos of helicity 1/2 by the oscillating scalar field in general is not suppressed by any 
powers of Mp^, and therefore their production can be very efficient. The magnitude of the 
effect is model-dependent. For example, this effect does not occur in the simplest model of 
a single chiral multiplet with a quadratic effective potential m?^(fp /2. Meanwhile this effect 
is very strong in the theory with the effective potential A0^/ 4. Gravitinos in this theory are 
produced very quickly, within about ten oscillations of the infiaton field, with occupation 
numbers rtk ~ 1/2 over a large range of momenta k < y/Xcj), which is very different from 
usual perturbative production. 

This result may have important cosmological implications, since it may allow to 
rule out certain classes of cosmological theories. The nature of this effect resembles the 
well known fact that the longitudinal components of massive vector bosons at high energy 
behave in the same way as the Goldstone boson that was eaten by the vector field 
A similar effect is known to exist in the theory of technicolor A more direct analogy 
is the transmutation of gravitational interactions of gravitinos with helicity 1/2 to weak 
interactions found by Fayet |jTD|]. In our case the effect is nonperturbative, and its adequate 
interpretation is achieved by finding solutions of the gravitino equations in a nontrivial 
self-consistent cosmological background. 

In order to study conformal properties of gravitinos we reformulated the standard 
N=l phenomenological supergravity in an SU(2, 2|l)-invariant way, which makes the confor- 
mal properties of the theory manifest and explains how the conformal symmetry is broken. 
Our formulation describes arbitrary number of chiral and vector multiplets and is flexible 
enough to allow investigation of regimes where the superpotential W vanishes. 

In application to the theory of gravitino production, we concentrate on the simplest 
models with one chiral multiplet $ and arbitrary superpotential. We present classical 
equations of motion and constraints for the transverse and longitudinal gravitino in the 
expanding Friedmann universe interacting with the moving infiaton field. We use the gauge 
where the goldstino is absent. Then we solve classical equations for gravitino. This solution 
confirms the generic prediction from the SU{2, 2|l)-symmetric theory that the longitudinal 
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gravitinos are not conformal. 

We represent equations describing gravitino components with helicities 3/2 and 1/2 
in a form analogous to the equations for the usual spin 1/2 fermions with time-dependent 
mass. This allows to reduce, to a certain extent, the problem of gravitino production to the 
problem of production of particles with spin 1/2 after preheating Q]. 

Finally, we estimate the number density of gravitinos produced by the oscillating 
scalar field in several inflationary models, and show that in some models the ratio 11-^/2/ s 
may substantially exceed the bound 1^1^/2/ s < O(10~^^). 

A detailed account of our investigation will be given in a separate publication |Tl 



Here we will only outline the main points of our study and present the most interesting 
results. 



2 Supergravity Lagrangian and conformal properties 
of gravitino 

Fundamental M-theory, which should encompass both supergravity and string theory, at 
present experiences rapid changes. One may still expect that the low-energy physics will 
be described by the N=l d=4 supergravity [0 and address the issues of the early universe 
cosmology in the context of the most general phenomenological N=l super gravity- Yang- 
Mills- matter theory [|T^. 

We are interested in conformal properties of supergravity fields, which include var- 
ious spin fields, in the conformally flat FRW metric describing the early universe: 



(2.1) 



In particular, we will be interested in conformal properties of gravitino. Supergravity is 
not a conformally invariant theory, despite the fact that it has a long history of being 



derived using the superconformal tensor calculus |jTj] as a technical tool. Therefore it is 
difficult even to address this issue as the supergravity fields do not have specific conformal 
weights. To solve this problem, the idea is to view the supergravity theory as a gauge- 
fixed version of the conformally invariant theory describing the most general = 1 gauge 
theory super conformally coupled to supergravity [[l^]. The derivation SU{2,2\1) invariant 
Lagrangian of = 1 supergravity coupled to + 1 chiral multiplets (with complex scalars 
Xj and fermions Qj) and Yang-Mills vector multiplets (with gauginos A° and vectors VT") 
superconformally will be presented in 
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It has no dimensional parameters. It consists 
of 3 parts, depending, respectively, on a real function J\f, a holomorphic function W and a 
gauge group 2-tensor fap- Each of them is conformally invariant by itself. 



C 



superconf 



mx,x)]n + mx)]p + luixyxixfi 



(2.2) 



The statement of the 5'?7(2,2|1) symmetry of the action (^) includes, among others, the 
symmetries under the following set of local dilatations, with parameter a{x), for the metric 
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and gravitino, for the scalars and spinors of the chiral multiplets, and for the vectors and 
spinors of the gauge multiplet, respectively: 

A"' = et"("U". (2.3) 

The function of scalars Af{X, X) codifies the information on Kahler manifold. The holo- 
morphic function of scalars W(X) codifies the superpotential. They transform as follows 
under local dilatations: 

AfiX', X') = e^'^^^WiX, X) , W(X') = e^'^^^^WiX) , ^X') = ^X) . (2.4) 

The important term in the conformal action which allows us to distinguish between confor- 
mal properties of helicity ±3/2 and ±1/2 gravitino is the following: 

mne-' = ^^{X, X) [r + i,^R^ + e'^d^^e^ ■ 'y^)) +■■■■ (2.5) 



The gauge fixing of the local dilatation of the conformally invariant action presented in 
leads to the standard^ Poincare supergravity theory. 

The dimensionful constants, Plank mass and F- and D-terms, appear via the gauge 
fixed value of the conformal compensator superfield, which we call conformon. The original 
n + 1 complex variables Xj are split into one complex scalar conformon field Y, and n 
physical complex scalars Zi, which are hermitian coordinates for parametrizing the Kahler 
manifold in the Poincare theory. One defines 

Xj = Yxi{zi), (2.6) 

and the local dilatation takes the form in which only the conformon Y transforms and the 
physical scalars do not transform under the local dilatation 

F' = e"(^'V, zl = Zi. (2.7) 

In these variables the gauge-fixing of the dilatational invariance ( p.3|) , ( p.7|) is given by 

M{Y, F, z, z) = -3|F|2 exp (-\lC{z, z)) = -3Mj,, (2.8) 



where the first equation defines }C{z,z), which is the Kahler potential, and the second is 
the gauge fixing. Here Mp = Mpianck/ VStt ~ 2 x 10^^ GeV. Thus the conformon field Y 

^ In fact two different gauges can be used to fix the i?-part of the superconformal symmetry. With the 
first choice W = W* we get the N=l phenomenological supergravity depending on Q (with e^^ = 
Mp^e'^\W\'^), while the second one Y ~Y* (see notation below) gives the version which is non-singular in 
the limit of the vanishing superpotential W = 0. This version is closer to the one in |]l^ which was obtained 
by the superspace methods. The limit W = from the first version has been discussed at the end of |]l6| . 
However, we find the second version more suitable for cosmology. It will be presented below. 
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is frozen to \Y\ = Mpexp (^^}C{z , z)"^ , i.e. it becomes a functional of the Kahler potential 
and not an independent field. The theory becomes that of the Poincare supergravity theory 
with 

[Ar]De~' = ~^Mj, [r + tfj^R'^) +.... (2.9) 

Here i?^ = j'^P'^Vpip^, where Vp is a covariant derivative. The local dilatation of the metric 
and gravitino in ( p.3|) is not compensated anymore by the local dilatation of scalars, 



Y = Y' ^ e"(^')F , Af{X, X) = U{X\ X') e2"(^W(X, X) . (2.10) 
The same happens with the W part of the theory: 

w(x) = Y^Mp'^w{z) =^ y\;(x) = yv{x') ^ e"'^''^yv{x) . (2.11) 

We have chosen here to give W mass dimension 3. 

Thus, after the freezing of the conformon field some part of the transformation 
cannot be performed and therefore some parts of the phenomenological supergravity La- 
grangian are not invariant under dilatations. One can try to change the dilatational weight 
for these fields to compensate the appearance of powers of Mp. However, this does not help, 
since the terms with derivatives on the conformon field are absent after the gauge fixing. 
The gravitino field equation which follows from the superconformal action is 



_ ^i^^^Q^ In + 7 ■ In + . . . = . (2.12) 

In the FRW cosmological problems only time derivatives of the scalar fields are important, 
therefore in ^^^dy In a/— A/" only the term "^^d^ In \J —M is relevant. After gauge fixing the 
conformal symmetry will be broken for configurations for which either 

7-V'7^0 or VoT^O. (2.13) 

Only such terms will be sensitive to the absence of the terms In ^ —H due to gauge 
fixing when —M = 3Mp. The gravitino in the general theory with spontaneously broken 
supersymmetry will be massive. The states of a free massive spin 3/2 particle were studied 
by Auvil and Brehm in |]l8l (see also for the nice review). A free massive gravitino has 
7 ■ = 0. Helicity ±3/2 states are given by transverse space components of gravitino, 
ijjf. Helicity ±1/2 states are given by the time component of the gravitino field ipQ. In 
cases when gravitino interacts with gravity and other fields, we will find that j ■ ip 0. It 
will be a function of tpQ. Thus the consideration of superconformal symmetry lead us to a 
conclusion that helicity ±1/2 states of gravitino are not conformally coupled to the metric. 
When these states are absent, the ±3/2 helicity states are conformally coupled (up to the 
mass terms, as usual). Thus the conformal properties of gravitino are simple, as it is known 

•^In [0, where an attempt to study conformal properties of gravitino has been made, it was assumed 
that 7 • ?/; = i/)" = and conformal symmetry of gravitino was deduced in the context of pure supergravity. 
Without scalars, however, pure supergravity does not support a cosmological background. In the presence 
of matter the assumption that j ■ ip = ip'^ = is not valid and conformal symmetry is broken. 
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for scalars: if the action has an additional term ■^^(p'^R, the massless scalars are conformal. 
If this term is absent, the scalars are not conformal. Note that both these statements are 
derivable from the super conformal action. We will see the confirmation of this prediction 
in the solutions of the gravitino equations below. 



As we already explained, our formulation starting with super conformal action [|TT 
provides flexibility in the choice of the form of = 1 phenomenological supergravity. If we 
take W = W* gauge for i?-symmetry, we get the action |jl3|| depending on the combination 
of the Kahler potential and the superpotential, called Q. Here we use the Y = Y* gauge for 
i?-symmetry and present the form of the phenomenological Lagrangian in which the Kahler 
potential and the superpotential are not combined in one function. This allows to avoid 
problems which sometimes appear when the superpotential W vanishes. The action can be 
written as 



R + i,.R^ + CsG 



torsion 



+ (Re U^) \-\F^,F^''^ - -pX^] + ii(Im Up) fcF^'^^ - 8, (A^TsT^A^ 



2^K 



Mp"e 



-mW* + {V'W)g-^i>{VjW)\ - i(Re 



+ |(Re u)^,^'^ [f:^ + f:;) 



1; • 



-e^/\V'Vm)xrXj + |i(Re /)-^"^Pa/^,X^A^ - 2MUa'g^'\''XJ 
+ lMp'e'^/\Vm)g-'/UA%X 



a 

R 



--Ai^R ■ 1X^\1X^ + U)X:X3>^l\l + h.c] 

+ M|.(?/ {\e-h^^-'^''i,,^,ijpx'l.Xi - i^,x' VX^ 
+ Ml {r1; - '^g,'g/) xYXkXi 



AMp2 ((Re /«^)A"7m75A^) - ^,Mp' P.iK^lg-'^ f.SjXlX 



^R 



+ |(Re f)-"^^ (/:,x.A^ - fa,^r^') {fpsXA' - fpsjX'^' 



(2.14) 



The L and R denote left and right chirality, e.g. \l = |(1 + 75)A, while for the x, the 
chirality is indicated by the position of the index: Xi is left chiral, while x* is right chiral. 
For the scalars, is the complex conjugate of Zi. The Kahler metric is gi^ , which is used 
also for covariant derivatives and Kahler curvature 



-pjk 



R 



M 



grd.T 



ki 
'3^ m ■ 



(2.15) 



Extra i indices on quantities, e.g. /^^ denote derivatives, here the derivative of fa/s with 
respect to Zi. Other covariant derivatives and notations are (antisymmetrization [/xz/] with 
weight 1, metric signature (— , +, +, +), notation = and 75 = 170717273) 



+ ^[.7.] A° , F^, = 2d[,W:^ + W^WJf,^ 



Pi ' 



2 c t ^ pa 1 ^ ~ ^ 1 '^At 



5m - w^;^" , 



A 



B 



{dilC)dz' - {d'lC)dzi + ?>W^{rl - r, 



(2.16) 



where 6a is a symbol for the transformation under the gauge group for all fields. For the 
conformon field Y and for the rest of the scalars we have 



(2.17) 



where r and ^ are n+1 holomorphic functions for every symmetry, such that the quantities 
in (12.21) are invariant. That determines also 



z,z*) + 3r*(z' 



(2.18) 



See flTI] for details. 

The appearance of Mp = \Y\e~^^^ in various places in this Lagrangian shows that 
the conformal symmetry is broken. One can rescale the fields with Mp so that they have 
standard kinetic terms. For our purpose it will be convenient to replace the scalar field 
by chiral fermions Xi by and similar for the gravitino, ipp. 



3 Gravitino equations 

In general background metrics in the presence of complex scalar fields with non-vanishing 
VEV's, the starting equation for the gravitino has in the left hand side the kinetic part 
and a rather lengthy right hand side which will be given in |Tl| . Apart of varying gravitino 



mass m = Mp'^e^W, the right hand side contains a chiral connection (see ( p.l6| )) 



and various mixing terms like those in the 4th, 5th and 6th lines of the phenomenological 
Lagrangian ( |2.14| ). For a self-consistent setting of the problem, the gravitino equation 
should be supplemented by the equations for the fields mixing with gravitino, as well as 
by the equations determining the gravitational background and the evolution of the scalar 
fields. 

Let us make some simplifications. We consider the supergravity multiplet and a 
single chiral multiplet containing a complex scalar field z = with a superpotential W 
and a single chiral fermion x- This is a simple non-trivial extension which allows to study 
gravitino in the non-trivial FRW cosmological metric supported by the scalar field. A 
nice feature of this model is that the chiral fermion x can be gauged to zero so that the 
mixing between ?/^^ and x iii ( 2.14| ) is absent. We also can choose the non- vanishing VEV 



of the scalar field in the real direction. Re $ = Im $ = 0, so that A^ = 0. The field 
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(j) = \pi Re $ plays the role of the inflaton field.0 Then from ( p.l4| ) we can obtain the 
master equation for the gravitino field 

+ ^V'm = (^M - j , (3.1) 
where gravitino mass m = m{(f){ri)) is given by 



W 



m = e^l^ ^ . (3.2) 



Gravitino equation ( |3.1| ) is a curved spacetime generalization of the familiar gravitino equa- 
tion (j^ + mo)"?/^^ = in a flat metric, where mg is a constant gravitino mass. 

The generalization of the constraint equations d^il)^ = and 'y^ipfj, = reads as 

+ ImYi^^ = , (3.3) 
fm^T^^M + m'a-2707V, = -^G^^Y^ , G^, = R^, - \g^,R , (3.4) 

where ' stands for a conformal time derivative 9^. It is important that the covariant deriva- 
tive in these equations must include both the spin connection and the Christoffel symbols, 
otherwise equation V^-^^ = used for the derivation of these equations is not valid. 

The last equation will be especially important for us. Naively, one could expect 
that in the limit Mp oo, gravitinos should completely decouple from the background. 
However, this equation implies that this is not the case for the gravitinos with helicity 1/2. 
Indeed, from ( p.4|) one can find an algebraic relation between 'y^ipo and Y'^f- 

7V0 = ^7Vi • (3.5) 

Here A is a matrix which will play a crucial role in our description of the interaction of 
gravitino with the varying background fields. If p and p are the background energy- density 
and pressure, we have Gq = Mp'^p, G\ = —Mp'^p6l, and one can represent the matrix A 
as follows: 

p - Sm'^Ml 2m'a~^Ml 
p + ?,m'^Ml ' p + 3mml. 



^ = ^ , o„....2 + ^0—^;-^ = ^1 + 70^ . (3.6) 



Note that in the limit Mp 00 with fixed ^ in z = ^Mp-^, one has m = Mp'^W . 
If W does not blow up in this limit, this matrix A is given by 

A = - + 7o , (3.7) 

P P 

where ' stands for derivative dt, and the relation between physical and conformal times is 
given by dt = a{r])dr]. In the limit of flat case without moving scalars, A = —1. 

■^Typical time evolution of the homogeneous inflaton field starts with the regime of inflation when (f> 
slowly rolls down. One can construct a superpotential W which provides chaotic inflation for (p > ^^p- 
When (f){t) drops below ~ Mp, it begins to oscillate coherent oscillations around the minimum of its the 
effective potential V {((>). 
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For definiteness, we will consider the minimal Kahler potential K, = zz* 



. In 



the models where the energy-momentum tensor is determined by the energy of a classical 
scalar field and $ depends only on time we have 



p = \^\^ + V, 



P 



Ml 



(3.8) 



We will use the representation of gamma matrices where 70 = diag (i, i, — i, — i). Then in 



]6D the combination A = Ai+iA2 emerges. For a single chiral multiplet we obtain \ A\ = 1. 
(One can show that |y4| = 1 for the theories with one chiral multiplet even if the Kahler is 
not minimal.) Therefore A can be represented as0 



A = — exp (21/ dt jj,{r]) 
\ J —00 



(3.9) 



Using the Einstein equations, one obtains for n (for minimal Kahler potential, and real 
scalar field): 



H = Mp^e^/^VVW + W) + 3{HMp^^ 



nil 



nil = m 



H = da 



-1 



(3.10) 



The expression for becomes much simpler and its interpretation is more trans- 
parent if the amplitude of oscillations of the field $ is much smaller than Mp. In the limit 
$/Mp one has 

/i = d^d^W . (3.11) 

This coincides with the mass of both fields of the chiral multiplet (the scalar field and 
spin 1/2 fermion) in rigid supersymmetry. When supersymmetry is spontaneously broken, 
the chiral fermion, goldstino, is 'eaten' by gravitino which becomes massive and acquires 
helicity ±1/2 states in addition to helicity ±3/2 states of the massless gravitino. 

The matrix A does not become constant in the limit Mp —>■ 00. The phase ( |3.9|) 
rotates when the background scalar field oscillates. The amplitude and sign of A change 
two times within each oscillation. Consequently, the relation between j^ipo and Y'^i ^^Iso 
oscillates during the field oscillations. This means that the gravitino with helicity 1/2 (which 
is related to ipo) remains coupled to the changing background even in the limit Mp 00. 
In a sense, the gravitino with helicity 1/2 remembers its goldstino nature. This is the main 
reason why the gravitino production in this background in general is not suppressed by the 
gravitational coupling. The main dynamical quantity which is responsible for the gravitino 
production in this scenario will not be the small changing gravitino mass m{t), but the mass 
of the chiral multiplet fi, which is much larger than m. As we will see, this leads to efficient 

^Initial conditions at inflation at 77 — + — cxi correspond to p = — p, m' = and A = —1, which gives 
/i(— 00) = 0. Alternatively, we can start with inflaton oscillations at 77 = 0, which defines the phase up to 
some constant. The final results depend only on fi. 
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production of gravitinos in the models where the mass of the 'goldstino' nonadiabatically 
changes with time. 

We shall solve the master equation ( |3.1|) using the constraint equations in the form 



y fj.-^j and (|3.3| ). We use plane-wave ansatz ~ e for the space-dependent part. Then 
ipi can be decomposed^ [|l9l into its transverse part ipf, the trace ^i^nd the trace k ■ ip: 



= ri + - 2^^(k ■ 7) j + [i^h - 27.(k ■ 7) j k ■ ^ , (3.12) 

where ki = A:i/|k|, so that Yipf = k^ipj = 0. We will relate "^'^ipi with ipQ and with k ■ 
so that, after use of the field equations there are two degrees of freedom associated with 
the transverse part ipj , which correspond to helicity ±3/2; and two degrees of freedom 
associated with 7*'^j (or ipo) which correspond to helicity ±1/2. 

For the helicity ±3/2 states we have to derive the equation for ipj . We apply 
decomposition ( ^.12| ) to the master equation (O) for /i = i and obtain]] 



a' 



y^'^^ + ^7° ± maj V^f = . (3.13) 

In the limit of vanishing gravitino mass, the transverse part ipj is conformal with a weight 
-1-1/2. The transformation ipj = a~^^'^'^J reduces the equation for the transverse part to the 
free Dirac equation with a time-varying mass term ma. It is well known how to treat this 
type of equations (e.g. see [Q). The essential part of "iff is given by the time-dependent part 
of the eigenmode of the transversal component ?/t(^); which obeys second-order equation 

+ {k'^ + nj, - in'j?j yr = , (3.14) 

where the effective mass is Qt = m{ri)a{r]). 

The corresponding equation for gravitino with helicity 1/2 is more complicated. We 
have to find k'^ipi and 7^1- The equation for the components k^ipi can be obtained from the 
constraint equation (p.3|): 



ik- ip = i 7o + 17 ■ k — ma ] 7^1 • (3.15) 

\ a J 

Combining all terms together, we obtain the on-shell decomposition for the longitudinal 
part 

ilji = ipf + iki'j-k + ^{^h - 7ik ■ 7)(^7o ± ma) j 'y^ipj . (3.16) 



^We use now "0^ with i = 1, 2, 3 for the space components of V'/j , while for gamma matrices 7* are space 
components of flat 7", and similarly for the index. 

similar equation obtained in ||] has a different coefficient in the term 1^7^^ since they have omitted 
the Christoffel symbols in the covariant derivative. One can still use their equation if one replaces the curved 
space gravitino vector -0^ for which the equation was derived by the tangent space vector = SaV'/i- 
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Now we can derive an equation for j^ipj. From the zero component of (|3.1|) we have 



3 of 

^tVo + ( T^rna + ik ■ 7 ) 



(7Vi) - ^7o7Vi • 



This equation does not contain the time derivative of ipQ. Substituting ipQ from 
( p.l7|) , we get an equation for j^ipi 



(d^ + B-ik- 77oi) 7V* = . 



where 



and 



A = - expl 270 / dtij,{r]) 



" 3a' 2 ma ^ j. 



(3,17) 
3|) into 

(3.18) 

(3.19) 

(3.20) 



We can spht the spinors 7*-?/' in eigenvectors of 70, 7^4 = + 6*-, and d± = |(1 =F 170) 7 Vi- 
From the Majorana condition it follows that 9±{k)* = ^C6z^{—k), where C is the charge 
conjugation matrix. In a representation with diagonal 70 the components 6± correspond to 
the 7o-eigenvalues ±i. Acting on (|3.18|) with the hermitian conjugate operation gives us a 
second-order differential equation on the 6^. We choose for each k a spinor basis Mi,2(fc) for 
the two components of 6^, and two independent solutions of the second-order differential 
equations fi,2{k,ri). The general solution is given by 

^+ = E a^^{k)Uk,v)up{k), 

a,l3=l 



9. 



A* k--f 
|A|2 P 



(d, + B)e^ = -C~' E a*-^{-k)f:{-k, vHi-k) . (3.21) 



The last equality determines reality properties of the coefficients. Here we represented B as 
Bi + 70-B2 and defined B = Bi + ii?2, by analogy with the definitions for the matrix A. By 
the substitution fa{k,r]) = E{ri)yL{Ti), with E = (— A*)^/2exp {— dr]Re B) , equation for 
the functions faik^rf) is reduced to the final oscillator-like equation for the time-dependent 
mode function yLijf)'- 

y'i+[e + nl-in'L)yL = Q. (3.22) 
Here ^ ^ 

o^'^VLl = fJi — -iJsin2 j jjidt + -m ^1 — 3 cos 2 j fidt^ . (3.23) 

In the derivation of ( p.22| ) it was essential that A has the form (|3.9|) . 

Finally we give the expression for the energy density of the longitudinal mode 

3 



a^Tn° 



4k2 



#(-k)7i 



i7 ■ k I ma - —70 



ma 7o 

a , 



ma + — 7o 7^V'i(k) . 
a / 



(3.24) 

-1, ""i^ipj ~ From To° we can define the occupation 
number of gravitinos of energy ujk at given mode Tq" = / d^kukn^, where is expressed 
through a bilinear combination of mode functions yiiv)- 



In the fiat spacetime limit Ai 
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4 Gravitino production 



One could expect that the gravitino production may begin already at the stage of inflation, 
due to the breaking of conformal invariance. However, there is no massive particle produc- 
tion in de Sitter space (i.e. as long as one can neglect the motion of the scalar field during 
inflation). Indeed, expansion in de Sitter space is in a sense fictitious; one can always use 
coordinates in which it is collapsing or even static. An internal observer living in de Sitter 
space would not see any time- dependence of his surroundings caused by particle production; 
he will only notice that he is surrounded by particle excitations at the Hawking temperature 
H/2TI. 

Gravitino production may occur at the stage of inflation due to the (slow) motion of 
the scalar field, but the most interesting effects occur at the end of inflation, when the scalar 
field (f) rapidly rolls down toward the minimum of its effective potential V{(f)) and oscillates 
there. During this stage the vacuum fluctuations of the gravitino field are amplified, which 
corresponds to the gravitino production (in a squeezed state). 

Production of gravitinos with helicity 3/2 is described in terms of the mode function 
yriv)- This function obeys the equation ( |3.14| ) with Qt = 'mcL, which is suppressed by 
Mp^. Non-adiabaticity of the effective mass i^riv) results in the departure of yriv) from 
its positive frequency initial condition e^''^, which can be interpreted as particle production. 
The theory of this effect is completely analogous to the theory of production of usual 
fermions of spin 1/2 and mass m [Q. Indeed, Eq. ( 3.14 ) coincides with the basic equation 
which was used in 0] for the investigation of production of Dirac fermions during preheating. 

The description of production of gravitinos with helicity 1/2 is similar but somewhat 
more involved. The wave function of the helicity 1/2 gravitino is a product of the factor E{ri) 
and the function i/Liv)- The factor E{ri) does not depend on momenta and controls only 
the overall scaling of the solution. It is the function yiiv) that controls particle production 
which occurs because of the non-adiabatic variations of the effective mass parameter QlIt]). 
The function yL^v) obeys the equation ( p.22| ) with the effective mass i^L('7), which is given 
by the superposition ( 3.23| ) of all three mass scales in the problem: /i, H and m. 

In different models of the inflation, different terms of Ql will have different impact 
on the helicity 1/2 gravitino production. The strongest effect usually comes from the largest 
mass scale n, if it is varying with time. This makes the production of gravitinos of helicity 
1/2 especially important. 

To fully appreciate this fact, one should note that if instead of considering super- 
gravity one would consider SUSY with the same superpotential, then the goldstino x (which 
is eaten by the gravitino in supergravity) would have the mass 9$9$iy, which coincides with 
H in the limit of large Mp, see Eq. ( |3.11|) . As a result, Eq. ( p. 22 ) describing creation of 
gravitinos with helicity 1/2 at ^ Mp looks exactly as the equation describing creation 
of goldstinos in SUSY. That is why production of gravitinos with helicity 1/2 may be very 
efficient: in a certain sense it is not a gravitational effect. (On the other hand, the decay 
rate of gravitinos T ~ rn? /Mp is very small because it is suppressed by the gravitational 
coupling Mp"^.) 
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This does not mean that one can always neglect terms proportional to H and m as 
compared to /i, and that production of gravitinos with helicity 3/2 can always be neglected. 
In order to understand the general picture, we will consider several toy models where the 
effective potential at the end of inflation has simple shape such a.s V ^ 0". We will not 
discuss here the problem of finding superpotentials which lead to such potentials (and 
inflation) at > Mp |2y], because we are only interested in what happens after the end of 
inflation, which occurs at ~ Mp. 

First consider the superpotential W = ^m^^'^. At <^ Mp it leads to the simple 



2 



quadratic potential V = -i^4> ■ The parameter in this case coincides with the infiaton 
mass m^. In a realistic inflationary model one should take ~ 10^^ GeV, which is equal 



to 5 X 10~^Mp I^Oj. Hubble constant during the field oscillations is given by , where 
00 (t) here is the amplitude of the field oscillations, which decreases during the expansion of 
the universe. The gravitino mass is given by m = -j^- 

Thus, at the end of inflation in this model, which occurs at ~ Mp, all parameters 
determining the behavior of the gravitino wave function are of the same order, fi ~ ~ 
H m. However, later the amplitude of decreases as 0o ~ ^'m^t ~ Tn'> where N is the 
number of oscillations of the field after the end of inflation 0]. Thus already after a single 
oscillation there emerges a hierarchy of scales, /x ~ ^ H ^ m. 

Since = const, after the first oscillation the parameter /i becomes nearly con- 
stant, the parameters H and m become very small, and their contribution to the gravitino 
production becomes strongly suppressed. As a result, the dominant contribution to the 
gravitino production in this model occurs within the first oscillation of the scalar field after 
the end of inflation. Each of the parameters /i, H, and m at the end of inflation changes by 
0{m^) within the time 0(m^^). This means that (because of uncertainty relation) graviti- 
nos of both helicities will be produced, they will have physical momenta k = 0{m^), and 
their occupation numbers will be not much smaller than 0(1). This leads to the following 
conservative estimate of the number density of produced gravitinos: 713/2 ~ 10"^m^. 

Now let us assume for a moment that all energy of the oscillating field transfers to 
thermal energy ~ within one oscillation of the field 0. This produces gas with entropy 

/ ^2 j^,j2 \ 3/4 

density s ~ ~ ( ^ | . As a result, the ratio 71,3/2 to the entropy density becomes 

~ 10- f ^V'' ~ 10-° . (4.1) 
s \Mp) ^ ' 

This violates the bound < 10~^^ for the gravitino with m ~ 10^ GeV by about 4 orders 
of magnitude. Thus one may encounter the gravitino problem even if one neglects their 
thermal production. 

In this particular model one can overcome the gravitino problem if reheating and 
thermalization occurs sufficiently late. Indeed, during the post-inflationary expansion the 
number density of gravitinos decreases as a~^. The energy density of the oscillating massive 
scalar field p = m|0o(t)/2 also decreases as a~^. But the entropy produced at the moment 
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of reheating is proportional to p^''^, so it depends on the scale factor at the moment of 
reheating as a'"^^^ . If reheating occurs late enough (which is necessary anyway to avoid 
thermal production of gravitinos), the ratio ^^^^ ~ 10~^°a~^/^ becomes small, and the 
gravitino problem does not appear. 

But this simple resolution is not possible in some other models. As an example, 
consider the model with the superpotential W = V^$^/3, which at <^ Mp leads to 
the effective potential A0^/4. The oscillations of the scalar field near the minimum of 
this potential are described by elliptic cosine, (f){ri) = ^cn{\/X(f)o, ^). The frequency of 

oscillations is O.8472-\/A0o and initial amplitude 0o — Mp p|. 

The parameter /i for this model is given by p = \/2X(j). It rapidly changes in the 
interval between and v^2A0o- Initially it is of the same order as H and m, but then H 
and m rapidly decrease as compared to fi, and therefore the oscillations of fi remain the 
main source of the gravitino production. In this case production of gravitinos with helicity 
1/2 is much more efficient than that of helicity 3/2. 

The theory of production of gravitinos with helicity 1/2 in this model is similar to 
the theory of production of spin 1/2 fermions with mass v^2A0 by the coherently oscillating 
scalar field in the theory Xcp'^/A. This theory has been investigated in [Q. The result can be 
formulated as follows. Even though the expression for contains a small factor v^2A, one 
cannot use the perturbation expansion in A. This is because the frequency of the background 
field oscillations is also proportional to -\/A. Growth of fermionic modes ( p.22| ) occurs in 
the non-perturbative regime of parametric excitation. The modes get fully excited with 
occupation numbers ~ 1/2 within about ten oscillations of the field 0, and the width of 
the parametric excitation of fermions in momentum space is about \/A0o- This leads to the 
following estimate for the energy density of created gravitinos, 

PS/2 ~ (VA0o)' ~ AV(0o) , (4.2) 
and the number density of gravitinos 

ns/2 ~ X'/'V'/\<Po) . (4.3) 

Now let us suppose that at some later moment reheating occurs and the energy 
density V^(0o) becomes transferred to the energy density of a hot gas of relativistic particles 
with temperature T ~ V^^'^. Then the total entropy of such particles will be s ~ ~ V^^^, 
so that 

^ ~ A3/4 ~ 10-10 . (4.4) 
s 

This result violates the cosmological constraints on the abundance of gravitinos with mass 
~ 10^ GeV by 4 orders of magnitude. In this model the ratio does not depend on 
the time of thermalization, because both 71^/2 and V{(f)o)^/4 decrease as a~^. To avoid this 
problem one may, for example, change the shape of V{(j)) at small 0, making it quadratic. 

The situation in models with V{(f)) ~ 0" for n > 4 is even more dangerous because 
the energy density of the oscillating field in such models decreases faster than a~^, so 
the entropy of the particles produced during reheating is suppressed stronger than by the 



14 



factor of a ^. The later reheating occurs in such models, the greater will be the resulting 



ratio 

s 



The most dangerous situation occurs in the class of inflationary models where the 
effective potential V{(j)) does not have a minimum, but instead monotonically decreases 
and becomes fiat at ^ oo. Such models have been studied recently by many authors 



21| , 1^, which gave them many different names: deflation, kination, and quintessential 
inflation. Following [Q, we will call them non-oscillatory (NO) models, which reflects an 
unusual non-oscillatory behavior of the scalar field after inflation. This behavior implies 
that the standard mechanism of reheating does not work in NO models. Therefore until 
very recently it was assumed that in such models all particles are produced gravitationally, 
due to the breaking of conformal invariance A typical model of such type has the 

effective potential which is given by V{(j)) ~ A0^/4 at < 0, and then it rapidly vanishes 
as becomes positive, so that Vi^cj)) — > at — > cxd. 

We will leave apart the question whether it is easy to obtain realistic versions of NO 
models in supergravity. For us it is only important that in such models the parameters /i, 
H and m change by 0(V^(/)o) = 0{H) during the time when the field (j) rolls from 0o ~ Mp 
to 0. Just as in one of the examples considered above, this should lead to production of 
gravitinos with number density which can be estimated as ~ 10~^H^. This number is of the 
same order as the number density of all other conformally noninvariant particles produced 
by gravitational effects in the scenario of Ref. Thus, barring the subsequent dilution 



of gravitinos by some late-time entropy release, one has 7x3/2/3 = 0(1), which contradicts 
observational data by 14 orders of magnitude. 

This problem of NO models can be resolved if one assumes that the scalar field 
interacts with some other particles a with a sufficiently large coupling constant g. This 



leads to production of particles in the context of the instant preheating scenario This 
mechanism is much more efficient than the gravitational particle production studied in |pT|] , 
and the entropy s of produced particles becomes much greater than 0{H^). This leads to 
a strong suppression of ^3/2/^ j^. 

But do we really have the gravitino problem in all of these models? In our investiga- 
tion we studied only the models with one chiral multiplet. This is good enough to show that 
nonthermal gravitino production may indeed cause a serious problem, but much more work 
should be done in order to check whether the problem actually exists in realistic models 
with several different chiral and vector multiplets. 

First of all, one should write and solve a set of equations involving several multiplets. 
Even in the case of one multiplet it is extremely difficult, and the results which we obtained 
are very unexpected. The situation with many multiplets is even more involved. One 
possibility is to consider the limit Mp 00, since the most interesting effects should still 
exist in this limit. 

But this is not the only problem to be considered. In the toy models studied in 
this section with W = |m(^$^ and W = |V^$^, the superpotential W and the gravitino 
mass vanish in the minimum of the potential at = 0. Then after the end of oscillations 
supersymmetry is restored, superhiggs effect does not occur and instead of massive grav- 



is 



itinos we have ordinary chiral fermions. In order to study production of gravitino with 
hehcity 1/2 with nonvanishing mass m ~ 10^ GeV one must introduce additional terms in 
the superpotential, and make sure that these terms do not lead to a large vacuum energy 
density. 

Models with one chiral superfield which satisfy all of these requirements do exist. 
The simplest one is the Polonyi model with W = a((2 — ^/3)Mp + $). One can introduce 
various generalizations of this model. However, potentials in all models of this type that 
we were able to construct are much more complicated than the potentials of the toy models 
studied in this section. In particular, if one simply adds small terms a + /5$ + ... to the 
superpotentials ~ $^ or one typically finds that V becomes negative in the minimum 
of the potential, which sometimes becomes shifted to the direction 02? where $ = {(f)i + 
1^2)/ V^- This problem can be easily cured in realistic theories with many multiplets, which 
is another reason to study such models. 

It would be most important to verify, in the context of these models, validity of our 
observation that the probability of production of gravitinos of helicity 1/2 is not suppressed 
by the gravitational coupling. We have found, for example, that the ratio for the 
gravitinos with helicity 1/2 in the model Xcp^/A is suppressed by A^^/^. This suppression 
is still rather strong because the coupling constant A is extremely small in this model, 
A ~ 10"^'^. However, in such models as the hybrid inflation scenario all coupling constants 
typically are O(10~^) [^. If production of gravitinos in such models is suppressed only 



by powers of the coupling constants, one may need to take special precautions in order to 
avoid producing excessively large number of gravitinos during preheating. We will return 
to this question and present a more detailed description of the effects discussed above in a 



separate publication [ITT 
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